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Abstract: The paper presents a joint subcarrier-pair based resource allocation algorithm in 
order to improve the efficiency and fairness of cooperative multiuser orthogonal frequency 
division multiplexing (MU-OFDM) cognitive radio (CR) systems. A communication 
model where one source node communicates with one destination node assisted by one 
half-duplex decode-and-forward (DF) relay is considered in the paper. An interference-limited 
environment is considered, with the constraint of transmitted sum-power over all channels 
and aggregate average interference towards multiple primary users (PUs). The proposed 
resource allocation algorithm is capable of maximizing both the system transmission 
efficiency and fairness among secondary users (SUs). Besides, the proposed algorithm can 
also keep the interference introduced to the PU bands below a threshold. A proportional 
fairness constraint is used to assure that each SU can achieve a required data rate, with 
quality of service guarantees. Moreover, we extend the analysis to the scenario where each 
cooperative SU has no channel state information (CSI) about non-adjacent links. We 
analyzed the throughput and fairness tradeoff in CR system. A detailed analysis of the 
performance of the proposed algorithm is presented with the simulation results. 

Keywords: cognitive radio; cooperative communication; resource allocation; proportional 
fairness; spectrum sharing 
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1. Introduction 

Cognitive radio technology (CR) has been proposed as a relatively new concept for improving the 
overall utilization of spectrum bands. This promising technology can allow the unlicensed secondary 
users (SUs, also referred to as CR users or CRUs) to access those frequency bands which are not 
currently being used by licensed primary users (PUs) in a given geographical area [1,2]. Cooperative 
communication technology [3] allows network nodes with single antennas to use other network nodes' 
antennas to transmit data, which can generate a virtual multiple-input multiple-output (MIMO) system. 
Cooperative spectrum sensing is a viable sensing technique to enhance spectral utilization efficiency of 
secondary users while ensuring the quality of service (QoS) of primary users [4]. In a CR network, 
SUs are allowed to transmit over the frequency bands of PUs as long as the resulting aggregate 
interference is kept below a certain threshold. This threshold is known as interference temperature 
constraint or interference power constraint [1]. As SUs can design power and subcarrier allocation 
strategies subject to such interference power constraints, the interference introduced to PUs is 
effectively controlled. A great deal of resource allocation algorithms and interference control strategies 
has been proposed for spectrum- sharing CR networks. For example, the optimal power allocation 
strategies to maximize the transmitted data rate of the secondary user with an effective protection of 
the primary user were studied in [5,6] for spectrum- sharing CR networks. 

Orthogonal frequency division multiplexing (OFDM) is an attractive modulation scheme for users 
in a CR system due to its flexibility in allocating resources among SUs. Since both SUs and PUs may 
exist in side-by-side bands, yet have different access technologies, mutual interference is the limiting 
factor for the performance of both networks. Thus, using of the classical subcarrier allocation and 
power loading algorithms, such as uniform power but variable rate and water-filling algorithms 
maximizing the transmission capacity of an OFDM-based conventional wireless network may result in 
higher mutual interference in the PUs' band. There is only one group of users in such a wireless 
network, i.e., PUs, for a CR system. 

According to the latest literature on resource allocation in cooperative communication [7-17], the 
relay users in the system do not transmit their own data and merely help other non-relay users transmit 
data. In some wireless applications such as cellular networks, however, each user has its own data to 
transmit so that it should allocate its total constrained power and subcarriers properly in transmitting 
its own data and relaying other users' data [18,19]. Tourki [19] focused on efficiency issues by 
studying how to maximize the total transmitted data rate in non-orthogonal amplify-and-forward (AF) 
cooperative scheme, which ignores the fairness among the cooperative users. According to [20,21], 
equal power allocation (EPA) among subcarriers was proposed to separate the user selection from the 
power of subcarrier. With EPA, the EPA-PRG (proportional rate greedy) [22] algorithm is proposed to 
maximize the system throughput while keeping the fairness. However, cooperative transmission 
technology isn't applied in this algorithm. In [23], a linear water- filling scheme (LWF-PI) was 
proposed. This algorithm maximized the overall transmitted data rate of the CR system while keeping 
the interference introduced to the PU bands below a threshold. However, the fairness among users was 
ignored. Chandrashekar et al. [24] proposed an algorithm which is capable of maximizing the total 
transmitted data rate and achieving a high proportional fairness index. However, this algorithm cannot 
be applied to the CR network where we must adjust the interference introduced to the PU bands below 
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a threshold. Tan [25] proposed a joint subcarrier and power algorithm based on Blotto games. This 
algorithm can achieve a good trade-off performance between fairness and efficiency in OFDMA-based 
cognitive radio network (CRN), but it cannot obtain the effectiveness of multiuser diversity for the 
SUs without ability to generate a virtual MIMO system. 

A novel scheme was presented in [26] for the allocation of subcarriers, rates, and power in 
orthogonal frequency-division multiple-access (OFDMA) networks. The resource-allocation problem 
was solved by decomposing it into a hierarchy of sub-problems. A joint subcarrier and power 
allocation algorithm was presented in [27] for cooperative MU-OFDM CR systems. In [28], a survey 
of resource allocation and scheduling schemes in OFDMA wireless networks was presented. 
Nader et al. in [29] considered the practical case in which only partial CSI for the wireless channel 
between the secondary base station and SUs is available at the secondary base station. They formulated 
the resource allocation problem in the secondary network as an optimization problem in which the 
objective was to maximize the weighted sum rate of the secondary users. A novel sub-channel and 
transmission power allocation scheme was proposed in [30] for multi-cell OFDMA networks with CR 
functionality. Tianxiang et al. in [31] discussed optimization over the relay assignment, subcarrier 
allocation, per node power control, and heterogeneous quality-of-service (QoS) provisioning. 
Sabit et al. in [32] investigated the performance of an OFDM-based CR spectrum sharing 
communication system that assumed random allocation and absence of the PU channel occupation 
information. Hong Xu et al. in [33] formulated a unifying optimization framework based on Nash 
bargaining solutions to fairly and efficiently allocate resources between primary and secondary 
networks, in both decentralized and centralized settings. As the optimal resource allocation scheme 
was highly complex, G. B. et al. [34] proposed a low complexity suboptimal subcarrier and power 
allocation scheme. They also proposed a suboptimal subcarrier allocation scheme that can guarantee a 
certain level of fairness among CR users. Naeem et al. introduced in [35] a hybrid heuristic algorithm 
for the relay assignment and power allocation problem which is a non-convex mixed-integer non-linear 
optimization problem, and this problem is generally non-deterministic polynomial-time (NP)-hard. 

In this paper, a joint subcarrier-pair based resource allocation algorithm in order to improve both 
efficiency and fairness index is presented first. The definition of fairness is borrowed from the 
networking literature. In contrast with [36], where large channel fluctuations are intentionally created 
with "dumb" antennas for long-term proportional fairness resource allocation, this paper proposes a 
subcarrier-pair based resource allocation algorithm to maintain proportional rates among SUs for each 
channel realization, which ensures the rates of different SUs to be proportional in any time scale of 
interest. By formulating the resource allocation and pairing problem in this way, it will be shown that a 
high transmitted data rate for all SUs (even those with poor channel gains) can be achieved with low 
computational complexity. Moreover, we extend the analysis to the case in which each SU can only 
have access to CSI of its adjacent links. This is a more realistic scenario when network nodes are 
mobile and the timely CSI cannot be exchanged between cooperative users. Consequently, each user 
can only have access to statistical CSI of non-adjacent links. It is shown that the system performance 
deteriorates due to limited CSI but still outperforms that of equal power allocation scheme. The key 
contributions of this work are: 
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1. It is considered that SUs need to transmit their own data directly to the destination, and 
in the next phase they also help their partner forward the data received in previous 
phase to the destination. Simulation results show that in the same situations the system 
transmitted data rate by proposed algorithm is the highest than that by LWF-PI 
algorithm [23], EPA algorithm [21] and the Optimal Scheme [37]. 

2. The proposed subcarrier-pair based resource allocation algorithm ensures the rates of 
different SUs to be proportional in any time scale of interest, simulation results shown 
that a high transmitted data rate for all SUs (even those with poor channel gains) can 
be achieved. 

3. It is considered that SU has no CSI about non-adjacent link. In this case, we take full 
advantage of the statistical information of the non-adjacent links. 

Notation: In this paper, a variable with "underline" m denotes the temporary optimal values within 
each iteration process, the "double underline"«denotes the optimal value, and the "bar" i denotes the 
statistical average value. £"(B)denotes the expectation operator, and ■ denotes the optimal value when 
only partial CSI can be obtained by SUs. 

2. System Model and Problem Formulation 

We consider a hybrid network consisting of a primary network (PRN) and a cognitive radio 
network (CRN) as shown in Figure 1. The CRN consists of a CR access point (AP) and 2K SUs. The 
PRN and CRN co-exist within the same geographical area. The access mechanism/modulation format 
in SUs' band is OFDM. Our focus is mainly on the uplink radio resource allocation in the CRN. The 
SUs are trying to find the opportunity to access to the AP. 



Figure 1. A cooperative MU-OFDM CR uplink system. 




According to [37,38], we also consider that the frequency bands of bandwidth B\, B 2 , B L which 
have been occupied by L PUs are sensed by the CR system and known to SU transmitters. Every two 
SUs form a cooperative partner and they are relay node for each other. As shown in Figure 2, the kth. 
(1 <k< K) cooperative partner consists of two SU transmitters, kl and k2. As is assumed in [23,37,38], 
we consider the same side-by-side CR radio access model. The unoccupied bandwidth sensed by SUs 
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for opportunistic spectrum access is located on each side of L PU bands as shown in Figure 3. The 
available bandwidth for CR transmission is divided into N subcarriers based on OFDM system. It is 
considered that the access mechanism/modulation format in PUs' band is not known to the CR system 
and the bandwidth for each CR subcarriers is AfHz. Some symbols are shown in Table 1. 



Figure 2. Model for cooperative transmission. 

Node I 




Figure 3. spectrum access model of cognitive radio system. 
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In general, there are three instantaneous fading gains in the uplink transmission scenario shown 
in Figure 1 : 

(1) The gains between the SU's transmitter and SU's receiver or AP for the nth 
subcarrier denoted as h£t , , respectively. 

(2) The gains between the SU's transmitter and Zth PU's receiver, denoted as h^'^ . 

(3) The gains between the Zth PU's transmitter and the SU's receiver or AP, denoted as 
h**,h%, respectively. 

The channel gains are modeled as independent zero-mean complex Gaussian random variables, 
where ki denotes the Zth SU in Mi cooperation partner and pi denotes the Zth PU band. According 
to [39], it is considered that these instantaneous fading gains are perfectly known at the SU's 
transmitter. Specifically, we assume that the SU's receiver can estimate channel gains h^'^j and 

/ipj^ and report to the CR transmitter. In Section 3.2, we will study the case where the instantaneous 
fading gains of the non-adjacent links are not perfectly known at the SU transmitter but the statistical 
information of the non-adjacent links are known at the SU transmitter. Moreover, it is assumed that 
primary receiver can estimate the channel h^'^ which is reported to the SU transmitter through a 

common control channel. 
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Table 1. Table of Symbols. 



Symbol Definition 



K Number of cooperative partners in the CR network 



N Number of subcarriers 



Af Bandwidth of a subcarrier 



Ts Length of a slot 



Pkii " Transmission power of kl th SU on subcarrier n in the first frame 

Pki2 2( " } Transmission power of kl th SU on subcarrier n in the second frame 

Transmission power of k2 th SU on subcarrier n in the first frame 
Pk2™' n) Transmission power of k2 th SU on subcarrier n in the second frame 

, .„,„ r , -s the channel gain of the communication link from the kith SU to the kj th SU 

user on the nth subcarrier 



, ss ,„ the channel gain of the communication link 

ki,0 

from the kith SU to AP on the nth subcarrier 



ISP." 
"ki.pl 



the channel gain of the interference link 

from Zth PU to kith SU user receiver on the 72th subcarrier 



^ ps , n the channel gain of the interference link 

from Zth PU to kith SU user receiver on the nth subcarrier 



^ ps ,„ the channel gain of the interference link 

from /th PU to AP receiver on the nth subcarrier 



i z k\ ' z k2' z o >' ~ 1>2,3,4} me ac jditive noises at the corresponding node 

{ $1 ' $.'2 ' $o ) ' 1 = 1' 2, 3, 4} the interference introduced by the PUs into corresponding node 



2.1. Cooperative Transmission among SUs 

The scenario of a three-node DF diversity model is considered, where one source communicates 
with one destination assisted by one half-duplex relay, as shown in Figure 2. One transmission period 
is divided into two consecutive frames. Communication takes place in two phases (listening phase Tl 
and relaying phase T2, the definition is according to the working state of relay user) for each frame. 
The power allocation scheme for kth cooperative partner on subcarrier n is shown in Table 2. The 
source node broadcasts its signal to relay and AP in Tl, whereas the relay and AP listen. The relay 
decodes the signal and forwards it to AP in T2. It is denoted that the subcarrier n in Tl is pairing with 
subcarrier SPl(n) in T2 for first frame, and pairing with subcarrier SP2(n) for second frame. In the first 
frame, klth SU receives data in this time slot while Hth SU transmits a symbol x kl (t) with 
power level P£f m on nth subcarrier in Tl. The symbol is received by node 0 (AP) and overheard by 

fc2th SU as: 

>Vo (0 = K;; 4 p m Vn)x n (o + zf (o + € (0 

, (i) 
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Table 2. Power allocation scheme for kth cooperative partner on subcarrier n. 





Tl in First Frame 


T2 in First Frame 


Tl in Second Frame 


T2 in Second Frame 


su« 


p«,SPl(«) 
^1,1 


0 


0 


pn,SP2(n) 




0 


p",SPl(n) 
/ t2,l 


pn,SP2(n) 

1:2,2 


0 



During this interval, the kith. SU decodes its overheard signal as x kl (t) and transmits it to the AP on 
SPl(n) subcarrier in T2 with the power level P^fP (n) . Then the AP receives the signal as: 



WO = KT n) ^^^) + 4 2) (0 + €\t) (2) 



In the second frame, the roles of fclth SU and k2th SU are reversed. Similarly, kith SU transmits a 
symbol ^(f) with power level fJ. B 2 ' s 2 P2( '° on the nth subcarrier in Tl. The symbol is received by node 0 

(AP) and overheard by kith SU as: 

WO = VZ^S^xM + 4 3) (0 + &f\t) 

i (3) 

wo = K s 2U p nT n) ^) + 4f (o + C(o 

In T2 of second frame, the AP node receives the noisy signal which is relayed by kith SU with the 
power level P^ (n \i.e.\ 



WO =C 2( " ) >/^%(0+4 4) (0+C(0 (4) 

2.2. Mutual Interference between PU Bands and CR Users 

In the MU-OFDM CR system, due to the coexistence of PUs and SUs in side by side bands, it is 
necessary to consider the mutual interference between PUs and SUs. There are two types of 
interference in the system. One is introduced by the PUs into the SUs band, and the other is introduced 
by the SUs into the PUs' band. In what follows, we provide brief description and mathematical models 
for interference between SUs and PUs. 

2.2.1. The Interference Introduced into PUs by SUs 

CR interference is introduced into the PU spectrum by CR out-of-band (OOB) emissions. OOB 
emissions arise as a result of transmit pulse shaping such that a portion of the CR radiated power in a 
vacant subcarrier is leaked into neighboring bands occupied by the PUs. According to [23], the 
interference factor which is the integration of the power density spectrum of the nth subcarrier across 
the Ith PU band, and can be written as: 

W<U =l K P i 1 Ts \ , I ( -^r ) d f (5) 

where Ts denotes the symbol duration, d„i denotes the distance in frequency between the nth 
subcarriers of SU band and Ith PU band, and 5/ represents occupied bandwidth by Ith PU. It can be 
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shown from Equation (5) that the interference to PU band is related to the distance between SU band 
and PU band. 

2.2.2. The Interference Introduced into SUs by PUs 

The interference introduced into kith SU and AP node transmitting in nth subcarrier by Ith PU can 
be denoted as Jpi ik i, Jpi f0 , respectively. According to [37], the interference value Jpi ki , Jpi f0 can be 

written as: 



J n p m(d nl ,P PU )=\hZ I 2 f/' a/ E{I N (w,0 PU (en)}dw 

1 2 

J n pl ,o( d m . P PU ) =l Kilo |2 E{1 « (W ' *™ > dw (6) 



2 

2 



, , , ,-wm 1 / , ™^sin(vv-<:W /2 S , 

= — -f <p PU (em . . ) 

2nM sm(w-£-)/2 

where w represents the frequency normalized to the sampling frequency, ^{/^ • )}is the power 
density spectrum of the PU signal after M-fast Fourier transform (FFT) processing, 0 PU (e JW ) is the 
power density spectrum of the PU signal. The PU signal has been taken to be an elliptically filtered 
white noise process with amplitude P PU . 

According to [40], using a relay is advantageous when: 

I h ss,n I 2 I U ss ' SFl ("> I 2 I U ss ' n I 2 

frame 1: min( , ) > 



2 V jn 2 V iSPl(n) 2 V jn 



1=1 1=1 1=1 

I iss,n i2 I 7.w,SP2(«) i2 i iss,n i2 

(7) 

f o w kl,kl \Kl,0 1 v. I ^fc2.0 1 

frame 2 : min( - L , ) > ^ 

^ 2 _lV T" ^ 2 _lV |-SP2(«) 2 V jn 

1=1 1=1 1=1 

in selective DF mode. It is considered that the link of source node— relay node and link of relay 
node— destination node are better than that of source node— destination node. Otherwise, the relay 
keeps idle on subcarrier n in the relaying phase for x" k2 (t) orx* kl (t) . In this paper, we just consider the 
case in which the relays keep working on each subcarrier, i.e., the Equation (7) is always true. 

3. Optimization Problem Formulation 

In this section, we analyze the joint optimization of subcarrier-pair based resource allocation 
algorithm for OFDM-DF based on full CSI and partial CSI, respectively. We are interested in how 
each SU allocates its power properly across its own data and its relayed data so as to maximize the 
system transmitted data rate while maintaining reasonable fairness between SUs. The optimization 
problem is formulated firstly and then solved in the dual domain. It is assumed that the PUs have a 
constant-rate, constant-power transmission, while the SUs are capable to adjust transmit power over 
different fading states based on the CSI of the CR network. We study a type of constraint imposed 
over the secondary transmission to protect the PUs by limiting the interference introduced to the PUs 
below a threshold. 
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3.1. Resource Allocation and Subcarrier Pairing Scheme Based on the OFDM-DF 

The CR AP combines the received signals from the source node in Tl and the relay node in T2 
through the maximal ratio combining. The transmit power is adjusted in each SU's transmitter. 
According to [27] and [41], when the link of source node->relay node transmission is successful 
for entire DF process, the transmission rate of kith SU and &2th SU at n subcarrier in relaying 
mode, which is connected via the Shannon capacity formula, can be shown as i^ pl(n > and /^ SP2(n) ; 
respectively: 

. , pfi.SPl(n) pn,SPl(n) pn,SPl(n) 

r k r n) =^iog 2 {i +I mn(ifc- i 2 — »± ,uq; P — ^ + I 2 — ^ )} 



p»,SPl(») 
r k\,\ 


»■ n kl,0 


pn.SPl(n) 
|2 ~*1,1 










(=1 


pn,SP2(n) 
,(2,2 


»' n k2,0 


pn,SP2(n) 
|2 **2,2 


L 








1=1 



1=1 



7 ™ = ^log 2 {l + min(lft- I 2 — ,1ft- I 2 — "4 + 1^™ I 2 ^ )} 



(8) 



where cr 2 denotes the Additive White Gaussian Noise(AWGN) variance. Here, it is assumed that all the 
channel gains are constant during two frames and the link between cooperative partners are symmetric, 
i.e., h s k l'l 2 = ft^'^for all k. The factor 1/4 in Equation (8) results from the fact that the transmission takes 
four slots in the cooperative scheme. 

Let P£ff 1 ^ + P^i 1 ^ = P™: SPl ( n \ i = 1,2 where i denotes the ith frame. This formula means 

that the average of the transmit power of the source node P^ P ' W and that of the relay node P^f'^ is 
constrained to be p£- SPl( - n \ which is the allocated power on subcarrier n at the source node for direct 

transmission. According to [42], the solution to this problem is the transmitted data rate and it is 
maximized when: 

p«,SPl(n) p«,SPl(n) p«,SPl(n) 
I ni,» |2 "ft 1,1 _i iss.n |2 £*1J , i r .s.s.SPl(n) |2 "k2,\ 

i n h 1 1, 0 i , i tTj, i n i j ' i n, ~ n i . 



kl,k2 ' L kl,0 ' L ' k2,0 ' L 

^ 2 +E^2 ^ 2 +E^,o ^ 2 +E4 p r 

/=i ;=i ;=i 

pn,SP2(n) pn,SP2(n) pn,SP2(n) 

I I 2 =1 /Co I 2 + 1 C P2W I 2 — — L ( } 

^ 2 +E^,« ^ 2 +E^,o - 2 +ZC > 

;=i /=i /=i 
for all \<n<N,l<k<K 



Let: 



I fcf'" I 2 I h ss,n I 2 I /z'"-" I 2 I A" " I 2 

1 ,l k\,k2 1 n 1 "H,0 1 n 1 "*2,(> 1 n 1 "t2,*l 1 



X(l,i2 — L 'Tkl.O ~ L 'Tk2,0 L '7k2,k 

^+11^2 ° 2+ Hj;,,o ^ 

i=i 1=1 i=i 1=1 

Together with P£f m + /f 2 'f (n) = P k n f m , i = 1, 2 we can obtain that: 

y SPK») y » _ r » 

pn,SPl(n) _ / £2,0 p«,SPl(n) p«,SPl(n) _ / kl,k2 I klfi pn,SP\(n) 

r kl,\ ~ n _ n SVKn) r k\ '^2,1 ~ n _ n . SPl(n) k\ 
fk\,kl /kl,0~ r fk2,0 lk\,kl I k\S> / k2,0 

y" -y" v SP2 <") 
pn,SP2(n) _ / k2M I k2,0 pn,SP2(n) pn,SP2(«) _ / klfl pn,SP2(n) 

"k\,2 ~ n _ n SP2(n) k2 > *k2,2 ~ „ _ „ SP2(n) ~kl 

7k2,k\ I i2,0 / *1,0 Ik2,k\ 7k2,0 + /k\fi 



(10) 



(11) 
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Denote ^ spl(n) , ^" 2 SP2(n) as the equivalent channel gain given by: 

n SPl(n) « V SP2(«) 

n,SPl(n) _ ' k\,k2/ £2,0 n,SP2(n) _ / k2,kU klfl 

"« ~ „» _„« , SPl(n)'"*2 - „ _ „ SP2(„) (12) 

/ kl,k2 fklfi 'k2,0 I k2,k\ f k2,0 f klfi 

By now, we can unify the transmitted data rate as: 

t (13) 

rn,SP2(n)/p«,SP2(n)N _ 1 M n,SP2(«) p n,SP2(«) i 

i it2 \ r k2 >~ ^ lu fe2l 1 " 1 " 'Afc2 J 

The proportional fairness is used as the optimized objective to develop subcarrier-pair based 
resource allocation algorithm in order to maintain a balance between system efficiency and fairness. 
According to [43,44] and the theorem in [45], we can readily deduce that there exists one unique 
proportionally fair allocation which can be attained by maximizing the objective function 

K N 

*F = 'YJ2l (0 kPkn^ki ¥l(n) + 7^" 2 SP2< " ) ) over the feasible set, where co k is the weighting factor to make the K 

k=l n=l 

cooperative partners achieve the desirable transmitted data rate. Besides, we should keep the 
instantaneous interference introduced to the PUs below a certain threshold. The constraints include the 
aspects of satisfying the maximum power and interference constraints as well as the minimum rate 
requirements. Therefore, the resource allocation problem can be formulated mathematically as given 
in Equation (14). Constraint CI corresponds to the subcarrier allocation constraint that each subcarrier 
n only can be allocated to one cooperative partner. C2 and C3 define that the sum of all the 
transmission powers of a particular SU on different subcarriers can't be greater than the maximum 
allowed limit for that particular SU. C4 ensures the cumulative interference from all SUs and through 
all subcarriers on a particular PU should not be greater than the interference limit set. C5 ensure that 
each SU can obtain the minimum rate requirements. This constraint precludes the possibility of 
multiple SUs simultaneously transmitting at the same subcarrier: 

max v=±±a> kPkn {itr n) +i:r {n) ) 

k=l n=l 
K 



subject to CI: ^ Pkn <l, Pkn G{0,l},\fn 



C2: ± Pkn p k r m 4' C3: tp^r in) * 4 

«=i 2 n=l 2 (14) 

C4: it/UfSTX, +Pk n 2T n) s n k2, Pl +PkT (n XZ n) +rg H ' ) sgS ) te iaP 

k=l n=l 

c5: ±pji" k r n) +i:r (n) }^R k 

n=l 

for all / = {1, 2,..., L},£ = {1, 2,..., /:} 

where K denotes the number of cooperative partner, k denotes the Mi cooperative partner, L denotes 
the number of PUs, / denotes the Ith PU, ptn is the allocation indicator that equals 1 when the nth 
subcarrier is allocated to the kth cooperative partner and 0 otherwise, P, is the transmit power budget 
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for each cooperative partner, the factor of 1/2 for the terms of PJ2 which results from the fact that it is 
a normalization for the transmissions within the duration of a frame, Ith (r> denotes the maximum 
allowable interference level at the Ith PU receiver, Rk is the minimum transmitted data rate for kth 
cooperative partners. 

The optimal solution to Equation (14) can be found by performing an exhaustive search with 
computational complexity 0(K N Z) [46], where K N is the number of possible subcarrier allocations and 
Z is the complexity of a power allocation algorithm for each subcarrier allocation. To reduce the 
exponential computational complexity, a suboptimum resource allocation algorithm with less 
computational complexity is developed in the following. The dual decomposition approach is used to 
solve the problem. The dual problem of Equation (14) can be formulated as: 

K N N K K D N 



min max* = Il^/U^f™ +/?f ^^^-^P^^-^P^f 



3t (0) >3t a) > p,p k=ln=1 n=x fc=1 k=1 2 



n=l 



^(2)^(3)^(4) 



k=l 1 n=\ 1=1 k=\n=\ 



(15) 



K N 

p«,SP2(«) r.SP2(«) p «,SPl(«) r,SPl(n) 11 ,V 3< 4 >rV/! ( J «. spl (n) , r«,SP2(«) s p -, 
+r k\,2 ^kl,pl + r k2,l ^k2,pl ll + Zj A k i-ZjPkn^kl + 1 k2 > K k* 

k=l n=\ 

subjectto /gf « >0,1 = {CX\^\^\^}>0 
n s {1,2,. .., N},{i,j}s {1,2}, ks {1,2..., K} 

where i and j denote the ith and j'th SU of the kth cooperative partner, respectively. The values of 
{A®\Aj c 1) ,A£ l) ,Af i \A}, 4) } are the introduced Lagrange multipliers. In the future using, we denote a 
vector shown in Equation (16): 

1 (0)_ M (0) ,(0) 3 (o)-,t ,(d _ r3 (i) id) ji'if -rl<'> 3 (1 h r 

A, — y/l^ , A2 , — , sl N J , A, — [/Lj ,Aj , — , sl K J , A, — , — , A K J 

The Equation (15) can be decomposed into two layers of sub-problems. In the lower layer, we can 
get K sub -problems: 

N 



( l m X^, -k {2) X y> X 4) ) = max V L f plM ' SP2M p kn 

n=l 



subjectto Pkn e{0,l},Vne{l,2,...,N},Vke{l,2,...,K} 

n - I n \ p _ I pn.SPl(n) p«,SP2(«) pn,SPl(«) p«,SP2(«) \ > () 
V ~ yPkn I ' r ~ \ r kl,l ' r kl,2 > r k2,l ' r k2,2 J- u 

where (17) 

L 

rn,SPl(n),SP2(«) / p«,SPl(n) p«,SP2(«K _ / rn,SPl(n) rn,SP2(«K _ V" 1 y (3) r p«,SPl(n) o n p«,SP2(«)on 
L /c V^l >^2 + l k2 > Zj / '■'Nfcl.l +/ A:2,2 3 /c2,/>/ 



1=1 



pn,SP2(n) r.SP2(n) p «,SPl(«) r.SPl(n) -. ; (0) ; (1) pn,SPl(«) ; (2) p«,SP2(«) , ; (4) / rn,SPl(«) , rn,SP2(«) s. 
+ I kl,2 ^kl.pl +r k2,l ^k2,pl i~ A „ ~ A k ^kl ~ A k r k2 +A k \h\ +1 k2 ) 

V n = 1,2,- ••, N,\fk = 1,2,- -,K 
We suppose thatF^is the maximum value of the objective function in the lower layer. The master 
problem in the upper layer could be expressed as: 
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,<W>^ ^ 1 M ^ (18) 

subject to {^ 0) ,^,Af,^ 3 >,^ 4) }>0 

Since a dual function is always optimized by first optimizing some variables and then optimizing 
the remaining ones. We define a subcarrier pairing parameter J3 n m e {0,1} that takes 1 if the nth 
subcarrier in Tl is pairing to mth subcarrier in T2 and 0 otherwise. We first optimize the primal 
variables with the assumption that dual variables {A,^ 0) ,A^ 1} ,A k 2) ,A^ (3) ,A ( k 4) } are given. The resource 
allocation and subcarrier pairing process can be divided into four stages: 

(a) Allocating the optimal power factor {P k n {f 1(n \P k n {f 2(n \P k ^ Kn) ,P k n ^ n) } for SUs. 
P i "'j P1(n) and P k n 22 1{n) imply the power used for self -data transmission, respectively. 
P k {f 2(n) and P k n 2\ l{n) imply the power used for partner-data transmission, respectively. 

4. Allocating the optimal set of subcarriers Q.k for Mi cooperation partner, i.e., obtaining 
the optimal subcarrier allocation factor p kn . 

5. Optimal pairing process for the subcarriers which are allocated to Q.k, i.e., allocating 
the optimal subcarrier pairing factor fi n , m . 

6. After the temporary optimal primal variables have been obtained in each 
iteration process, we would find the temporary optimal dual variables 
{X^ 0) , , A} 2) , A^ 3) , A^ 4> } , which can minimize the objective function H(P) as shown in 

Equation (18). 

3.1.1. Power Allocation Algorithm 

^ tI ^ i n)^n)^ I ^(n) fI i^n)> ) = I n^m +I ^n) ^ if we make p ^ = 1, the power allocation 
can be determined in a water-filling fashion. Taking derivatives of £,^ SP1 <") SP2 <") with respect to 

p«,SPl(n) pn,SP2(n) 

C)T,i,SPl(«),SP2( n ) ( pn ,SPl(n) p „,SP2(nK cjpn,SPl(n),SP2(n) L 

U H \ r kl >_£k2 > _ (m , l(4)N^fH ;(D V 3(3)r/»«.SPl(n)c" _l_ F".SP1(") oSPl(n) n _ r> 

SP1(n) -{Ct) k +/l k ) SP1(n) \ 2-1*1 L/1 *i.i ^ki, P i +p k2,i ^ki.ni J~ u 

0/ M 1=1 

^rn,SPl(n),SP2(n) ( pn,SPl(i>) pn,SP2(n)\ ^pn,SPl( n ),SP2(n) L 

\<°k +A k ) P p n ,SP2(„) A k ~ 2-i A l \- U kl,2 k\,pl + tS k2.2 ^k2,pli~ U 



^jn,SP2(n) 

where 



dp nm dP k2 



SPlfn) SP2(«) 
,n,SPl(n) _ / k2,0 c n,SP2(n) _ /H,0 



(19) 



kU ~ n _ n SPl(n) ' *2,2 " „ _ „ SP2(„) 

/ kl,k2 /kl,0^/k2,0 I k2,k\ I k2,0 T / kl,0 

j-jn,SP2(n) _ Yk2,kl ~ ^2,0 z^n.SPHn) _ ^Hj2 ~ ^1,0 

«,2 - . „ +y SP2(«)' *24 ~ » _ y » +y SPl (n) 

/ fc2,fel /*2,0 T / H,0 /Jtl,t2 / /tl,0 T / /c2,0 

Taking derivatives of /?;- SP1 w. SP2 <"> = /£ SP1 <"> + 1"*™ with respect to P k f pl(n) , P k n > SV2{n) , and according 
to Equation (19), we can get: 
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gjjB,SPl(!>),SP2(B) 

~, p n,SPl(«) 
^ll 



A/ 



B,SP1(I>) 

lk\ 



41n21 + / 7 ;; SP1( " ) J P; i ' SP1( " > 



^«,SP1(«) 

«,+4 (4) 



g^n,SPl(n),SP2(n) 



A/ 



B,SP2(B) 

'Ik2 



/^n,SP2(n) 
^k.2 



(20) 



4 In 2 1 + V n k 2 nW P k n 2 m <»k + A <4) 



Q n,SP2(n) _ 1(2) .V 0 (3) r n «,SP2(n) C SP2(«) H n,SP2(n) c « 
&k,2 ~\ + Zj A I t U n,2 * kl.pl + tS k2,2 ^k2,pl\ 



Together with the constraint p^ spi ^ n \ pn,sP2(n) > q ^ tem p 0rar y optimal solution can be obtained: 



pn,SPl(n) _ 
r kl 



p«,SP2(«) _ 
"k2 



0, 



0, 



A/ 
41n2 



7*i 



n,SPl(n), 3 (4K ra n,SPl(n) 



Cjn,SPl(n) n,SPl(n) 
A/" n,SP2(n) r , o(4K c,n,SP2(n) 



n,SP2(n)„B,SP2(n) 



t,2 %2 



(21) 



The temporary SU's transmit power can be obtained for the given dual variables: 

p«,SPl(n) _ jn,SPl(n) pK,SPl(«) p«,SPl(n) _ pn,SPl(«) p«,SPl(n) 

11.1 — A jfcl,l £jj > 1:2,1 — **2,1 1l 

p«,SP2(K) _ r»n,SP2(«) pn,SP2(n) p«,SP2(n) _ p«,SP2(n) p«,SP2(n) 

11.2 — ^1,2 r t2 » 12,2 — "*2,2 *_k2 



(22) 



3.1.2. Subcarrier Allocation Algorithm 

The subcarrier allocation constraint is that each subcarrier is allocated to no more than one SU 
cooperative partner, which prevents mutual interference among SUs. According to Section 3.1.1, we 
can get a temporary optimum power {P£f 1(n \P^ l(n \P£f 2(n \P k n ^ (n) } . We substitute this 

temporary optimum power vector into the objective function l;' sp1( " )SP2<m) (P t " spl(M) , P^ 2 ' SP2(n) ) and the 
objective function i?; spl(n)SP2(M) (P; i spl( ' 1) ,p; 2 SP2(n) ) to obtain the temporary max value Lf pl( " )SP2(n) and 
^«,spi(«),sp2(«) ^ reS p ec ti ve iy Taking account of resource fairness, we can formulate the optimization 

problem of subcarrier allocation as: 

max%) = XZA„C!!!! + ^ 

P k=\ n=l ^^(^j/ 7 ^«,SPl(n),SP2(«)^2 

M n=l " (23) 

K 

subject to p kn e {0, 1} , £ p kn < 1, Vk, n, 

k=\ 

where p = {p k J,l< k < K,l< n < N 

where the value C is the weighting factor to balance the total transmitted data rate and fairness index 
among SUs. The bigger the value C is, the greater the fairness can be obtained, and otherwise the 
greater the transmitted data rate is. It is a mixed binary integer programming problem that is difficult to 
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solve. To reduce the exponential computational complexity, a suboptimum subcarrier allocation 
algorithm with less computational complexity is developed. The pseudo-code of subcarrier allocation 
algorithm can be described as follows: 

(1) Initialization 

(i) Make ^=0,0, =0,Vk,n 

(2) For n = 1 to N 

(i) For k = 1 to K 

a) Make pu n = 1 

b) Assign V n k = Hp) according to Equation (23) 

c) Make p kn = 0. 

(ii) Assign k* = argmaxV^ 



(iii) Assign p . =1,£\- =^\« u { n }>A„ =0,Vfc ^k 

According to Section 3.1.2, we can obtain a temporary subcarrier allocation vector 
P = (An'l - k<K,\<n<N) for the given dual variables. 

3.1.3. Subcarrier Pairing Algorithm 

The pairing constraint is that each subcarrier m in listening phase only pairs with at most one 
subcarrier n in the relaying phase. We assume that the pairing for deferent frames is not the same. The 
pairing process of the subcarrier allocated to kth cooperation partner can be expressed as: 



Constraints C6, C8 and CIO correspond to the pairing constraint that each subcarrier n in listening 
phase only pairs with one subcarrier m in the relaying phase in the first frame. Constraint C7, C9 and CI 1 
correspond to the pairing constraint that each subcarrier n in listening phase only pairs with one 
subcarrier m' in the relaying phase for the second frame. We can obtain the temporary optimal m, m' 
for any n as: 



Frame 1 : max £ £ fi„ F ki m ; Frame 2 : max £ £ fi^ F k ? 




subject to C6: X /? nm <1,C7: £ P n m , <l,VneQ, 



(24) 



C8 : X /?„,„, < 1, C9 : £ < LVm e Q„ Vm' g Q, 



C10:/?„, me {0,l},Cll:/? n , m ,e{0,l} 



m = argmax {/;j m }, rri_ = argmax {/; 2 m ' } 



(25) 



m 



in 



That is: 



Frame 1 : fi n m = 1, SPl(n) = m, J3 n m =0 Vm ^ m 
Frame 2 : J3 n r ^ = 1, SP2(«) = n/, J3 n m , = 0 Vm V 



(26) 
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The subcarrier pairing scheme can be shown as following: 
For k = 1 to K 

(a) A = n k ,A' = n k , A, m =0,^=0,V«,m,m'eQ ( 

(b) While A^0,A'^0 

(i) fieA,n'eA' 

(ii) Find m e A, m' e A' satisfying /£ m > /£' , /,"' m ' > /;•'' ,V/ e A,V/'e A' 

(iii) Assign 

SPl(n) = m, J3 n m = 1, A = A - {n} - {m} 

SP2(n') = m\ J3 n , m , = 1, A' = A' - - {m'} 
Through Section 3.1.3, we can obtain a temporary optimal subcarrier pairing vector /? nrn , (3 nm ', 
1 < {n, m, m'} < iV for the given dual variables. 

3.1.4. Optimizing the Dual Variables 

The optimal values of dual variables can be achieved iteratively by the sub-gradient method 
as follows: 

K 

; (0),i+l _ , ; (0),i (i) r, V H+ 

A <lw+1 ={A (iy -^[^-Z^nT SP1( ' !> ]} + 

^- «=1 

={2f-«-<4-t^p;™]} + (27) 

K=l 

/r iV 

j(3),i+l_M(3),i „(i)rf f L(i) Wo /p«.SPl(n) c« pn,SP2(«)r.n pn,SP2(n) C SP2(«) p n,SPl(n) r.SPl(n) yi^ 

A l -{A, -a 3l [lth ~ ^Z^Pkn^m s n, P i + ^2,2 \ 2 ,„/ +/ h,2 Vp( + / *2,i \ 2 , P / All 

*:=1 «=1 

a< 4) "' +i = { 4 4W - < (t a» (^ SP1W + nr n) ) - ^ i r 



where{«Q„,«j ! fc ,«2A' a 3/ ' a «}^ s a smai l positive step size for the z'th iteration, z is the iteration number 

A 

and(x) + =max(x,0) , with the appropriate step sizes, the iterations are convergent. The remaining issue 
is how to determine the step size {a^, a j l) k , a ( 2 l) k , a ( 3 l ] ,a^'l} . Clearly, performing a line search at each 
iteration process perform well. For a given current iteration X' ={/l^ 0)! ,/^ (1) ' , ,/^ <2) ' ! ,/l / (3) ' I ,/l it (4) ''} and a 
search direction d (,) , we compute step sizea*' (,) by: 

a 0,n' a ],k' a 2,k' a 3,l ' a 4,k 

yfc=l n=l £=1 

yt=l /=1 fc=l 
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Within each iteration process, the power allocation vectors can be updated respectively by 
Equation (22), the subcarrier allocation vectors can be updated respectively by subcarrier allocation 
algorithm shown in Section 3.1.2, the subcarrier pairing vectors can be updated respectively 
by subcarrier pairing algorithm which is shown in Section 3.1.3, with the updated value 
{/l ^ (0) '^/l ^ <1) •^>^f ) •^/l ^ (3) •^2f ) • , }. Therefore, the dual variable X 1 will converge to the dual optimum las 

i — >• oo and the temporary primal optimum variable will also converge to the primal optimum value 
after several iterations, i.e., temporary vectorf/f/f 1(B) ,/gf X{n \P k ;f 1{n) ,P k n ^ in) } will converge to optimal 

YGCtor{P k ^ Kn \P^ w ,P^ m ,P^ (n) }, temporary optimal vector p = [ Pkn ,l < k < K,\ < n < N} will 

converge to the optimal vector p = {p kn , 1 < k < K, 1 < n < N) , temporary vector P nm ,/3 n m , will converge 

to the optimal vector P nm ,/3 n ^ . 



3.2. Resource Allocation with Partial CSI 



If full CSI can be achieved at the Hth and kith (1 <k<K) SU transmitter, the optimal subcarrier 
allocation vector t , subcarrier pairing vector (3 nm (3 nmi , the power allocation vector 



{P k n ;f Kn \P k \f 2(n \P k fP {n \P k n z s ^ {n) } can be determined simply by the proposed algorithm which is 

shown in Section 3.1. However, the practical case in which only partial CSI of the wireless channel 
between the secondary base station and SUs is available have to be considered. The CSI of 
non-adjacent link may be undesirable and even unavailable when the SUs are mobile. In this section, 
we investigate the optimal resource allocation in SU cooperation network with partial CSI at each 
transmitter. Specifically, the kith. SU transmitter has full CSI of its adjacent links h^'^jfi j), 

but only statistical CSI of non-adjacent link/z/^f . We assume that the SUs can know h^' ki , , i.e., 

hpi'^hpiQ are still available at SU transmitter, and the link between cooperative users is symmetric for 

simplicity, i.e., K£ =Km • 

Under these assumptions and according to [47], the objective function: 

v=ii;^A.{^ ) (^" ) )+^ ) (^r ( " ) )} (29) 

k=\ n=\ 

can be rewritten as: 



Y = 22^A.{^r* ) (^ ) ) + ta/S ,HC " ) (flr o,) )} (30) 



where I n k f P[<n) and /;f 2(n) are defined as: 



tb,SP1(b) , p n,SPl(n) \ _ T? r tk,SP1(«) , p n,SPl(«) ^, 
l k\ )- t hss.nl 1 k\ ^ F k\ )J 

n k2,0 



r n,SP2(n) , n n,SP2(n)\ r- r T«,SP2(n) , n «,SP2(ii) »i 

hi y F i ) = t h ss.nU k 2 y^ki )i 



(31) 



In order to seek the optimal power allocation solution, we derive the explicit expressions for 



^.SPKn) and/ ».SP2 (n ) described as . 
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Let/z^' 0 " =( ) v x ki0 ,i = l,2, where d kj0 denotes the distance between the Mm SU and the AP, v is 



l ki,0 



the path-loss exponent, xia,o, (i = 1,2) is the normalized complex Gaussian random variable distributed 
as CN (0, 1), then at high signal-to-noise ratio (SNR): 



W 1 



(1- 



Af p".SPK») 1 W+W pn,5P\(n) n W,W 

y^,SPl(n) _ A/ r ffcl /H,t2"l ^ V M q Al,<r2 VK l ^ VK 1 x 

fcl /llO n n n n " " 

4mZ 7h,£2 ~~ 7*1,0 7fcU2 ~Ykl,0 Ykl,k2 ~Tkl,0 Yk\,k2 ~7kl,0 

(32) 



-)} 



Yk2,k\ Yk2,0 



4in2 r; 2 , H -r; 2 , 0 

where Wi,W2 in Equation (32) are defined as: 



(1- 



w 2 



Ykl.kl Yklfl Yk2,k\ Yk2 



-)} 



0 



(33) 



Proof: It is sufficient to show /^ spl(n) ; and 7 t " 2 SP2( " ) can be derived in exactly the same way, then 



r n,SPl(«) 



l k\ 



p r jn,SPl(n) i _ p r A f ln „ n , „i,SPl(n) p «,SPl(«) 
jt2,0 "*2,0 4 



tl )1 



n SPl(n) 
-F r , n , Ykl,klYk2fi 

" liTo 4 §2 y» v" + r SP1W ' H 

* 2 ' () ^ ^U2 _ ^l,0 + ^2,0 



pn,SPl(n)^ 



F 1 

n k2,0 



Yk\,k2 



I , m,SP1(«) ,2 
1 n k2fl 1 



Af 



a 2 + 



Z-, J pi,o 



log 2 (l + - 



/=i 



Y'k\k2~Yk\,0 + l 
/=1 



■ , M ,SP1(») ,2 *1 
1 n k2,0 1 



p n,SPl(nh 



Af , 

4 



*2,() 



„,« rv— 2 , \' r« \j2v i-l i „ i2 D «,SPl(nh 
Ykl,k2K<7 +2_ l J pl,<y> d k2fl} 1**2,0 1 *H ) 



log 2 (l + - 



(^l,fe2 ~ Yklfl + ~ 



I X (t2,0 ' 



C^ 2 _l V rSPl(n)s ,2v 
( +L J plfl > d k2,0 



(34) 



The third equality holds since we assume each transmission block is long enough to undergo 
different channel realizations as argued in [48]. Next, let: 

W { =[(a 2 + XC ) )<or 1 J=l %0 P (35 ) 



i=i 



then according to [49], we can obtain: 
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,n,SPl(n)_ F r ,n,SPl(n)-._ 4f F 

4 "fc2,0 



A/ f+oo. 



log 2 



1 + 



v" W vP"' SP1(n) 



y*U2-y*i,o+Wiy 



1 , ^U2^3<^ 
y*U2-yti,o+Wiy_ 

g -yf"^"V t y 2 Ty l(y " U2 -^ 1|0) ^ 



V 
Jo 



A/ f "MO 



41n2 



f 

Jo 



y + 



w, 



Yk\,k2~Yk\fi 



-}dy 



(36) 



y+- 



Af <-e w > 



41n2 

-A/ 



7k\k?~7k\(\ „» „» 7k\,k2-Yk\fi 
" khk2 " kX '° 7k\fi-7kU2 y 

f Wi — dy + e 

J-CO y 



7kl,0~7kl,k2 



Pkf^YkikiW^ rP k f mn) 7kkiW^ e l_ 

J-CO -y 



7kl,k2 7kl,0 



-e 



Yk\fl-YkU2 e ^ m 7 n kim^ 
41n2 ' W X "'^"V^i+W, 



P ., Ykl,0~Ykl,k2 



where £V(0) denotes the exponential integral function defined as = f e x x 'dx and £X0 can be 

J GO 

expanded asymptotically as: 

J -ATT' & e 1 1 

M(0 * - (1 + X ii T^) - - [1 + ~ + o(~)] (37) 
t r-f x t t t 

Then using the assumption that SNR is high, we readily obtain Equation (32). This completes 
the proof. 

We can use the resource allocation and subcarrier pairing algorithm which is proposed in 
Section 3.1 to solve the corresponding optimization Equation (14) once again. When partial CSI can be 
achieved by SUs, the optimal subcarrier allocation vector p , subcarrier pairing vector ^ nrn , /? n rn ' 

and optimal power allocation vector P£f P1 ^ , p^ SP2 ^ can b e obtained. According to Equation (11), we 
can get 

Yk2,Q 



*ki,i *h%" 



>n,SPl(n) | 



-L 



+ oo 



,o( v n _ v n . SPUn) kl 
\Ykl,k2 Ykl.Q ^ Yk2,0 



pn,SPl{n) - y , 

n _ n , w v r kl e u y 

o Yki,k2 Yki.o + vv iy 



Ykl,k2 Ykl.O 

Yki, k2 - Kfci.o + w iy. 



e y dy 



/fel.O /fel,fe2 



, ■ s /- + oo /• 

_ p n,SPl(n) I I . 



e 1 - 

— dt 
t 



(38) 



fc2 



,n,SPl(n) 



,,n _ v n fcl 
Kfci,/c2 Kfci,o 
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Then we also can get: 




p n,SP2(ri) p n,SP2(n) 



1 \ D n,SPl(ri) 



p n,SP2(ri) 
' "kl,2 



Yk2,kl Yk2,0 



w 2 



y n,SP2(ri) 
k2 



(39) 



Intuitively, the system performance would be degraded due to limited CSI though this scheme 
which does not require each cooperative node to have full CSI of nonadjacent link, as verified in the 
numerical simulation. 

3.3. Comparison with Classical Resource Allocation Algorithms 

Several existing schemes are compared with the proposed subcarrier-pair based resource allocation 
algorithm in terms of the system transmitted rate and fairness, respectively. These existing schemes 
include the following: 

LWF-PI-without-SP: The subcarrier is allocated according to the channel gain. The messages 
transmitted on subcarrier m at the source node will be retransmitted on the same subcarrier m at the 
relay node. The power is allocated according to LWF-PI algorithm [23] on each subcarrier. 

EPA-without-SP: The subcarrier is allocated according to the channel gain. The messages 
transmitted on subcarrier m at the source node will be retransmitted on the same subcarrier m at the 
relay node. The power is allocated equally on each subcarrier. 

LWF-PI-with-SP: The subcarrier is allocated according to the channel gain. The messages 
transmitted on subcarrier m at the source node will be retransmitted on subcarrier n, which is selected 
by proposed subcarrier pairing algorithm, at the relay node. The power is allocated according to 
LWF-PI algorithm [23] on each subcarrier. 

EPA-with-SP: The subcarrier is allocated according to the channel gain. The messages transmitted 
on subcarrier m at the source node will be retransmitted on subcarrier n, which is selected by proposed 
subcarrier pairing algorithm at the relay node. The power is allocated equally on each subcarrier. 

Optimal-Scheme-with-SP: The subcarrier is allocated according to the channel gain. The messages 
transmitted on subcarrier m at the source node will be retransmitted on subcarrier n, which is selected 
by proposed subcarrier pairing algorithm at the relay node. The power is allocated according to 
Optimal Scheme [37] on each subcarrier. 

Partial CSI: the difference between the previous five algorithms and Partial CSI scheme is that the 
SU has only imperfect CSI of non-adjacent link under partial CSI scheme. In this case, the objective 
function is formulated by statistical methods. 

4. Simulation Results 

We have studied asymmetric or linear network with all SUs of kth cooperative partner located in the 
same line. Specifically, Hth SU and the destination AP are fixed at (0, 0), (1, 0) respectively, and £2th 
SU is located at id, 0), 0 < d < 1, without loss of generality. The results for path loss exponent v = 2 are 
presented and all channels are modeled as Rayleigh flat fading with AWGN. Some simulation 
parameters are shown in Table 3. 
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Table 3. Simulation Parameters. 



Parameter 


Value 


Number of cooperative partners K 


4 


Number of PUs L 


2 


Number of subcarriers N 


20 


Length of a slot Ts 


4u s 


value of amplitude PP V 


10 x 10~ 3 W 


Ith il) (/= 1,2) 


2.7 x 10~ 3 W 


average channel power gain 


10 dB 


A/ 


0.315 MHz 


B 2 


1 MHz 


B 2 


2 MHz 



Simulation results are presented in this section to verify the performance of the proposed 
subcarrier-pair based resource allocation algorithm. In our simulations, the CVX, a package for specifying 
and solving convex programs, is used to solve formulated optimization resource allocation problems. 

The channel gains ^^,/?^' 0 ",^J p ",^^",/z^Q used in this paper are assumed to be Rayleigh fading, 

since the channel fading gains for different realizations of channel gain can be different, an average 
transmission capacity of 10,000 independent simulation runs is considered. And individual fairness 
index is defined as [46]: 

{ttp, K SPl(n) ' SP2(n) } 2 

fair = ^f^ (40) 

*Z(Za A"' sp1( " XSP2( " ) ) 2 

k=l n=l 

The fairness ranges between 0 and 1. The higher the value fairness is, the more fair the throughput 
distribution among SUs is. 

4.1. The System Transmitted Data Rate Obtained under Resource Allocation Algorithms 

As shown in Figure 4, the achievable maximum CR system transmitted data rate is plotted versus 
the power budget P t . The upper curve denotes the transmitted data rate by proposed subcarrier-pair 
based resource allocation algorithm with full CSI. It can be noted that the proposed subcarrier-pair 
based resource allocation algorithm achieves the highest transmitted data rate under a given total 
power constraint. The transmitted data rate achieved using the proposed algorithm is the highest 
among that using LWF-PI-with-SP, Optimal-Scheme-with-SP algorithm, LWF-PI-without-SP 
algorithm and EPA algorithm. The main reason is that the proposed algorithms can make full use of 
the entire available interference threshold, while the LWF-PI algorithms can only guarantee that the 
total interference is under the interference threshold as shown in Figure 5. We observe that only the 
proposed algorithm outperforms the compared existing resource allocation algorithms. Therefore, we 
can conclude that the proposed subcarrier-pair based resource allocation algorithm makes valuable 
contribution to system transmitted data rate. 
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Figure 5. Power Budget versus interference introduced to PU bands. 




4.2. The Fairness Index Obtained under Resource Allocation Algorithms 

Both the proposed resource allocation algorithm and EPA-with-SP algorithm exhibit the best 
fairness performance and LWF-PI-without-SP algorithm shows the least fairness as shown in 
Figure 6. Moreover, the fairness loss of EPA-with-SP allocation algorithm compared to the proposed 
subcarrier-pair based resource allocation algorithm is acceptable. In the LWF-PI-without-SP algorithm, 
LWF-PI-with-SP algorithm and Optimal-Scheme-with-SP algorithm, most power and subcarrier will 
be assigned to SUs with good channel conditions to improve system efficiency. Unlike the proposed 
subcarrier-pair based resource allocation algorithm and EPA schemes, these two schemes achieve 
significant performance improvement to ensure fairness among SUs. The fairness obtained by 
EPA-with-SP effectively approaches to 0.98. Due to multiuser diversity, the fairness attained by the 
proposed algorithm is above 0.95, while the achievable system transmitted data rate is higher than that 
of LWF-PI and EPA-without-SP algorithms. 



Sensors 2013, 13 



10327 



Figure 6. Power Budget versus Fairness index. 




0.65 1 ' ' 1 1 1 

0.005 0.01 0.015 0.02 0.025 

Power Budget ( Pt ) of CR user (in Watts) 



4.3. The Transmitted Data Rate of Each SU for the Resource Allocation Schemes 

As shown in Figure 7, the achievable transmitted data rate of each SU is plotted. The transmitted 
data rate of each SU, under the proposed subcarrier-pair based resource allocation algorithm, is the 
most balanced and stable among EPA, LWF-PI and Optimal Scheme [37]. Especially under the 
LWF-PI algorithm and Optimal Scheme [37], some SUs with bad channel conditions have lower 
transmitted data rates, which would result in unsuccessful communication with destination node. We 
can conclude that the proposed subcarrier-pair based resource allocation algorithm makes valuable 
contribution to balance transmitted data rates among SUs. 

Figure 7. Maximum transmitted data rate versus each SU (The power budget P, is 20 x 1CT 3 W). 




1 2 3 4 5 6 7 8 

the number of each SUs 



4.4. The System Transmitted Data Rate Obtained under Full CSI and Partial CSI 

As shown in Figure 8, the performance of proposed subcarrier-pair based resource allocation 
algorithm is presented for the case that each SU has no CSI of non-adjacent link. It is shown that the 
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system transmitted data rate decreases due to limited CSI compared to that of full CSI scenario. 
However, this scheme does not require the full CSI of non-adjacent link at each cooperative SU. 
Moreover, we show that the transmitted data rate of the proposed subcarrier-pair based resource 
allocation algorithm with partial CSI significantly outperforms that of EPA. The transmitted data rate 
is close to that of subcarrier-pair based resource allocation algorithm with full CSI. 

Figure 8. Power Budget versus sum transmitted data rate. 




0.005 0.01 0.015 0.02 0.025 

Power Budget ( Pt ) of CR users (in Watts) 



4.5. The Transmitted Data Rate Obtained by SUs under Different Distance 

Last but not least, we have also studied how the distance between cooperative SUs impacts on the 
system performance. Under the proposed subcarrier-pair based resource allocation algorithm with 
proportional fairness, Figure 9 shows that the system has comparatively better performance 
when d e [0.4,0.6] . This result provides a guideline for grouping and partner selection in user 
cooperative networks. 

Figure 9. Maximum transmitted data rate versus distance between two cooperative SUs. 




5 I i i i i i i i 
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5. Conclusions 

In this paper, we have developed a novel subcarrier-pair based resource allocation algorithm that 
maximizes the transmission data rate while the interference introduced to the PUs remains within a 
given limit. Using the proportional fairness as the optimized objective function, we can improve 
proportional fairness of resource allocation and achieve substantial transmitted data rate gains. The 
sum of power constraint for the source and relay nodes are considered. Moreover, we extended the 
analysis to the case that the CSI of nonadjacent link is not available at cooperative SUs, and found that 
even in this case the proposed schemes perform better than the classical schemes. 

Compared to the existing resource allocation algorithms which are introduced in the paper, our 
algorithm considers that each secondary relay user has its own data to be transmitted. Simulation 
results have shown that, either in improving the system throughput or in improving fairness of resource 
allocation, the proposed subcarrier-pair based resource allocation algorithm offers the best 
performance conditions among several existing compared resource allocation algorithms under various 
power budgets, while keeping the interference introduced to PU bands below a specified threshold. 
Besides, the transmitted data rate of proposed subcarrier-pair based resource allocation algorithm 
obviously outperforms that of EPA scheme when only partial CSI can be obtained by SUs. The 
contribution by taking full advantage of the statistical non-adjacent links channel information is 
demonstrated clearly in the simulation results. In addition, we notice that the system efficiency loss of 
partial CSI scheme compared to that of full CSI scheme is acceptable. 
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